Percolation with Multiple Giant Clusters 
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We study the evolution of percolation with freezing. Specifically, we consider cluster formation 
via two competing processes: irreversible aggregation and freezing. We find that when the freezing 
rate exceeds a certain threshold, the percolation transition is suppressed. Below this threshold, the 
system undergoes a series of percolation transitions with multiple giant clusters ("gels") formed. 
Giant clusters are not self- averaging as their total number and their sizes fluctuate from realization 
to realization. The size distribution Fk, of frozen clusters of size k, has a universal tail, Ft ~ fc -3 . 
We propose freezing as a practical mechanism for controlling the gel size. 
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Percolation was originally discovered in the context of 
polymerization and gelation Percolation has found 

numerous applications in physics 0, 0], geophysics ||, 
chemistry |fj , and biology |j| • It plays an important role 
in a vast array of natural and artificial processes ran ging 
from flow in porous media |8( and cloud formation [9|, Il0| 
to evolution of random graphs 0, 0] , combinatorial op- 
timization |l.'i| . algorithmic complexity amorphous 
computing, and DNA computing using self-assembly |l5| . 

We study the evolution of percolation using the frame- 
work of aggregation. An aggregation process typi- 
cally begins with a huge number of molecular units 
("monomers") that join irreversibly to form clusters 
("polymers"). At some time, a giant cluster ("gel") con- 
taining a finite fraction of the monomers in the system 
is born, and it grows to engulf the entire system. In this 
classic percolation picture only a single gel forms, but in 
many natural and artificial processes the system freezes 
into a non-trivial final state with multiple gels or even 
micro-gels ^3,0]. In this Letter, we show that aggrega- 
tion with freezing naturally lead to formation of multiple 
gels and that freezing is also a convenient mechanism for 
controlling the gel size. 

We analyze the simplest aggregation with freezing 
process where there are two types of clusters: active 
and frozen. Active clusters join by binary aggrega- 
tion into larger active clusters. The aggregation rate 
is prop ortional to the product of the two cluster sizes 
HHlEI; this is equivalent to the gelation model of 
Flory and Stockmayer where a chemical bond between 
two monomers joins their respective polymers [HH In 
parallel, active clusters may become frozen at a size- 
independent constant rate a. These frozen clusters are 
passive, that is they do not interact with other (active or 
passive) clusters. 

This process is conveniently studied using the rate 
equation approach. The density Cfc(t) of active clusters 
of mass k at time t (that is, made up from k monomers) 



satisfies 
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with m(t) the mass density of active clusters. The first 
two terms on the right-hand side describe how the cluster 
size distribution changes due to aggregation, and the last 
term accounts for loss due to freezing. The quantity m(t) 
is subtle. Generally, it equals the mass density of all 
active clusters including possibly giant clusters, but when 
there are no giant clusters, i.e., all clusters are finite in 
size, then m(t) = Mi(t) with the moments defined via 
M n (t) = (k n ) — J2k>i ^" c fc(*)- We are interested in the 
evolution starting with finite clusters only. 

The gelation transition. Initially, all clusters are finite 
in size, so m = Mi. The moments M n provide a useful 
probe of the dynamics. From the governing equation 
the second moment of the size distribution Mi obeys the 
closed equation dM2/dt = M 2 (M2 — a). For arbitrary 
initial condition, 
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There is a critical freezing rate a c = M2(0). For fast 
freezing, a > a Cl the second moment is always finite 
indicating that clusters remain finite at all times. In 
this case, there is no gelation. For slow freezing, a < a c , 
there is a finite time singularity indicating that an infinite 
cluster, the gel, emerges in a finite time |2fJ. The gelation 
time is 
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The gelation point marks two phases. Prior to the 
gelation point, the system contains only finite clusters 
that undergo cluster-cluster aggregation ("coagulation 
phase"). Past the gelation point, the gel grows via 
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cluster-gel aggregation ("gelation phase"). We analyze 
these two phases in order. 

Coagulation phase. Coagulation occurs for a > a c at all 
times or for a < a c when t < t g . From the mass 
density of active clusters satisfies dm/dt — —am, and 
thus, ordinary exponential decay occurs, 



m{t) = m(0)e 
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For concreteness, we consider the monodisperse initial 
conditions Cfc(O) = 5k,i- In this case M„(0) = 1 and 
consequently, a c = 1. The cluster size distribution is ob- 



tained using the transformed distribution, Cfc 
and the modified time variable 
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that increases monotonically with the physical time and 
reaches r — > 1/a as t — > oo. With these transformations, 
Eq. reduces to the no-freezing case 



i+j=k 

From the well-known solution of this equation |2lL [22 
the cluster-size distribution is 
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Generally, the size distribution decays exponentially at 
large sizes and the typical cluster size is finite. The gela- 
tion time © is simply r g = 1. Of course, no gelation 
occurs when a > 1 because r < 1/a < 1. Otherwise, 
as the gelation point is approached, t — * i g , the charac- 
teristic cluster size diverges A; ~ (t g — £)~ 2 . The gelation 
point is marked by an algebraic divergence of the size dis- 
tribution Cfc ~ (1 — a)fc~ 5 / 2 for large k. We note that the 
mass density decreases linearly with the modified time, 
m = 1 — ar, and that at the gelation point, the mass is 
simply m(r 9 ) = 1 — a. 

Gelation phase. Past the gelation transition, a giant clus- 
ter containing a finite fraction of the mass in the system 
forms. In addition to cluster-cluster aggregation, cluster- 
gel aggregation takes place with the giant cluster growing 
at the expense of finite clusters. In parallel, all clusters 
may undergo freezing and particularly, the gel itself may 
freeze. 

Formally, the size distribution Q generalizes to 
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with u(t) = f Q dt' m(t'). Statistical properties of the 
size distribution are derived from the generating func- 
tion c{z,t) — J2k>i kck(t)e kz that equals 



c(z,t) 
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FIG. 1: The total mass, the sol mass, and the gel mass versus 
modified time r for a — 1/2. These curves hold as long as 
the gel remains active. 



(6) where G ( z ) = J2k 



e kz is the "tree" function ]2c 



During the gelation phase, active clusters consist of 
finite clusters, the "sol", with mass s, and the gel with 
mass g. The total mass density of clusters, m = s + g, 
decays according to 



dm 
~~dt 
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The sol mass decays according to ds/dt = —gM 2 — as 
obtained from (JTJ using s = M\. The first two moments 
follow from the generating function, M\ — c(0) and 
M 2 = c'(0), and using the identity G'{z) = G/(l - G) 
|24|. yields M 2 = s/(l — sre at ). The evolution equation 
for the sol mass becomes explicit, 
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Equations (f 1 Cjf> and (|11[) are subject to the initial condi- 
tions m{tg) = s(t g ) = 1 — a. Once the masses are found, 
the formal solution |JSJ becomes explicit. Results of nu- 
merical integration of Eqs. (|10|l and (|ll(l are shown on 
Fig. H In the vicinity of the gelation transition, the gel 
mass grows linearly with time: 



g(t) ~ 2(1- a) 2 (t-t g ) 



(12) 



as 1 1 t g . The quadratic dependence on the freezing rate 
implies that the emerging gel is very small when a | a c . 
Thus, micro-gels, that may be practically indistinguish- 
able from large clusters, emerge. Moreover, the maximal 
gel size must be smaller than 1 — a. We conclude that 
freezing can be used to control the gel size, as gels of ar- 
bitrarily small size can be produced using freezing rates 
just below criticality. 

Multiple giant clusters. At any time during the gelation 
phase, the gel itself may freeze. This freezing process 
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FIG. 2: The mass density m versus time r. Shown are results 
of a Monte Carlo simulation with system size N = 10 s and 
freezing rate a — 0.1. The system alternates between the co- 
agulation phase and the gelation phase. In the former phase 
the mass decreases linearly according to Q such that deple- 
tion occurs at time r = 1/a. In the latter phase, the active 
mass decreases slower than linear according to 1101 and 1111 . 
The gelation phase ends when the gel freezes. 

is random: the gel lifetime T is a random variable that 
is exponentially distributed, P(T) — ae~ aT . Until the 
gel freezes, the system evolves deterministically, so the 
mass of the frozen gel is g{t g +T). When the gel freezes, 
the total active mass m(t) is discontinuous: it exhibits 
a downward jump (Fig. |2J). Given that the duration of 
the gelation phase is governed by a random process, the 
mass of the frozen gel is also random. It fluctuates from 
realization to realization, i.e., it is not a self-averaging 
quantity. 

When the gel freezes, the system re-enters the coagula- 
tion phase because all remaining clusters are finite. The 
initial conditions are dictated by the duration of the pre- 
ceding gelation phase and are therefore also stochastic. 
Nevertheless, once the initial state is set, the evolution in 
the coagulation phase is deterministic. The gel is frozen 
so it no longer affects the evolution and only cluster- 
cluster aggregation occurs. Let us assume that the gel 
freezes at time tf. Reseting time to zero, the first and 
the second moments are simply given by Eqs. J3J and 
©, respectively, with M n (0) replaced by M n (tf). Note 
that M n (tf) contains contributions from finite clusters 
only. A second gelation occurs if the freezing rate is suffi- 
ciently small, a < M2(t/). Otherwise, the system forever 
remains in the coagulation phase. 

Cyclic dynamics. The general picture is now clear: the 
process starts and ends in coagulation and throughout 
the evolution, the system alternates between coagulation 
and gelation. Once the initial conditions are set, the be- 
havior throughout the coagulation phase and throughout 
the gelation phase are both deterministic. Each gelation 
phase ends with freezing of the active gel. Since the dura- 




FIG. 3: The size distribution of the frozen clusters for a = 
1/2. The simulation results represent an average over 10 2 
independent realizations in a system of size N = 10 6 . 

tion of the gelation phase is random, the size of the giant 
clusters, and their number are both random variables. 
Generically, the system exhibits a series of percolation 
transitions, each producing a frozen gel, so that overall, 
multiple gels are produced. The random freezing process 
governs the number of percolation transitions as well as 
the size of the frozen gels. 

The magnitude of the second moment when the gel 
freezes determines whether a successive gelation occurs. 
Because the second moment diverges at the gelation 
point, there is a time window past the gelation time 
where the second moment exceeds the freezing rate. If 
the gel freezes during this window, another percolation 
transition is bound to occur. Therefore, the maximal 
number of frozen gels is unbounded. 

Monte Carlo simulations confirm this picture (Fig. EJ. 
In the simulations, we keep track the total aggrega- 
tion rate R a = N(M± — Ma)/2 and the total freezing 
rate Rf = <xNMq. Aggregation occurs with probability 
Ra/(Ra + Rf)i an d freezing occurs with the complemen- 
tary probability. A cluster is chosen for aggregation with 
probability proportional to its size. Time is augmented 
by At = l/(R a + Rf) after each aggregation or freezing 
event. 

Frozen clusters. We now turn to the frozen clusters. 
Their density is found from dFk/dt = ack- When a > a c , 
the density of active clusters over the entire time range is 
given by Eq. J7J, and therefore Fk(t) is found by simple 
integration. In particular, the final density is 

Fk ^= pTfci 7(M/«) (13) 

where 7(71, a;) = J Q X dyy 71 " 1 e~ y is the incomplete gamma 
function. At large sizes, the size distribution of frozen 
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clusters decays according to 
F k (oo) 
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-B(a)k] 



a — 1 
a > 1 



(14) 



where A = (27r) _1 / 2 a 2 /(o: - 1) and B = a' 1 + In a - 1. 

Similarly, we may obtain the density of frozen clus- 
ters formed in the first coagulation case. Integrating (JZJ) 
gives Fk{t g ) = l{k,k). This size distribution decays 



algebraically, Fk(t g 



However, enumerating the 



frozen clusters that are born past the first gelation point 
is more complicated since the deterministic patches are 
punctuated by random jumps. In fact, even the total 
mass of frozen clusters becomes a random quantity be- 
cause it is a complementary quantity to the final mass 
density of frozen gels. Interestingly, numerical simula- 
tions suggest that the tail behavior 
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is universal (Fig. OJ. Even the prefactor D appears to 
depend only weakly upon the freezing rate a. This shows 
that the final distribution of frozen clusters is dominated 
by the coagulation phase. Indeed, in each such phase, 
the tail behavior is fc -3 , and furthermore, in the gelation 
phase, large clusters are more likely to be consumed by 
the gel. These results suggest that freezing leads to an 
additional non-trivial critical exponent 7 = 3. 

In conclusion, we have found that freezing leads to mul- 
tiple percolation transitions. The system evolves in a 
cyclic fashion alternating between coagulation and gela- 
tion. Depending on the freezing rate, the system may 
form no gels, a single gel, or multiple gels. Most statis- 
tical characteristics are non-self-averaging because they 
are controlled by the random freezing of gels. 

Freezing provides a practical mechanism for controlling 
gelation. It may be used to engineer micro-gels of desired 
size by implementing variable cooling rates. Slow freezing 
followed by rapid freezing can be used to produce gels 
of prescribed size, while near-critical freezing produces 
micro-gels of arbitrarily small size. 

There are a number of interesting potential generaliza- 
tions of the present work. The most natural is percolation 
in finite dimensions. Percolation is analytically tractable 
only in two dimensions, though even in that situation it 
is not clear how to derive the exponent characterizing 
the tail of the final distribution of frozen clusters. One 
may also consider situations with different freezing mech- 
anisms |25| . particularly for finite and infinite clusters. 
We also solved for the case where infinite clusters freeze 
immediately [2(|. In this case, there is no breakdown 
of self-averaging and the final density of frozen clusters 
mimics the critical behavior of active clusters, 7 = 5/2. 
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